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Abstract 

Given a univariate polynomial, its abscissa is the maximum real part of its roots. 
The abscissa arises naturally when controlling linear differential equations. As a 
function of the polynomial coefficients, the abscissa is Holder continuous, and not 
locally Lipschitz in general, which is a source of numerical difficulties for designing 
and optimizing control laws. In this paper we propose simple approximations of the 
abscissa given by polynomials of fixed degree, and hence controlled complexity. Our 
approximations are computed by a hierarchy of finite-dimensional convex semidefi¬ 
nite programming problems. When their degree tends to infinity, the polynomial 
approximations converge in norm to the abcissa, either from above or from below. 
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1 Introduction 


Given a univariate polynomial, its abscissa is the maximum real part of its roots. When 
studying linear differential equations, the abscissa of the characteristic polynomial of 
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the equation is used as a measure of the decay or growth rate of the solution. In linear 
systems control, the abscissa function is typically parametrized by a small number of 
real parameters (the controller coefficients), and it should be minimized so as to ensure a 
sufficiently fast decay rate of closed-loop trajectories. 

As a function of the polynomial coefficients (expressed in some basis), the abscissa is 
a Holder continuous function (with exponent equal to the reciprocal of the polynomial 
degree), but it is not locally Lipschitz. As a consequence of this low regularity, numerical 
optimization of the polynomial abscissa is typically a challenge. 

For a recent survey on the abscissa function and its applications in systems control, see [B]. 
A detailed variational analysis of the abscissa was first carried out in [Sj. These ideas were 
exploited in a systems control setup in [3], using randomized techniques of non-convex 
non-smooth local optimization, however without rigourous convergence guarantees. 

In the space of controller parameters, the zero sublevel set of the abscissa function of the 
characteristic polynomial of a linear system is the so-called stabilizability region, and it is 
typically non-convex and non-smooth, see [8] where this set is approximated with simpler 
sets such as balls or ellipsoids. In [7j, ellipsoidal approximations of the stabilizability 
region were generalized to polynomial sublevel set approximations, obtained by replacing 
negativity of the abscissa function with positive definiteness of the Hermite matrix of the 
characteristic polynomial. 

This paper continues the research efforts of [8] and [7], in the sense that we would like 
to approximate the complicated geometry of the abscissa function (and its sublevel sets) 
with a simpler function, namely a low degree polynomial. The level of complexity of the 
approximation is the degree of the polynomial, to be fixed in advance. Moreover, we would 
like the quality of the approximation to improve when the degree increases, eventually 
converging (in some appropriate sense) to the original abscissa function when the degree 
tends to infinity. 

The outline of the paper is as follows. After introducing in Section 2 the abscissa function 
and some relevant notations, we address in Section 3 the problem of finding an upper 
approximation of the abscissa. In Section 4, we address the more difficult problem of 
approximating the abscissa from below, first by using elementary symmetric functions, 
and second by using the GauB-Lucas theorem, inspired by £5]. Explicit numerical examples 
illustrate our findings throughout the text. 


2 


2 Preliminaries 


Notation and definitions 

Let n G N and Q C M n be a compact semi-algebraic set on which a Bor el measure with 
support Q can be defined and whose moments are easy to compute. For simplicity, in this 
paper we choose Q = [—1,1]” = {q G M n : 1 — qf > 0,..., 1 — > 0}. 

Let ‘if(Q) denote the space of continuous functions on Q. Its topological dual is isometrically 
isomorphic to the vector space of signed Borel measures on Q. By Banach-Alaoglu’ 

s theorem BE!. the unit ball of ^C(Q) is compact (and sequentially compact) in the 
weak-star topology of 

Denote by M[g] rf the vector space of real polynomials in the variables q = (qi,..., q n ) of 
degree at most d. Let E[g] C M[g] be the convex cone of real polynomials that are sums of 
squares of polynomials and denote by S[g] 2c z its subcone of sum of squares polynomials of 
degree at most 2d. 


The abscissa function 

Consider the monic non-constant polynomial p <G M[s] defined by 

m 

P : s i )■ p(q,s) := 

k =0 

with s G C complex, q — {qi,...,q n ) £ Q and given polynomials pk € M[g] for k = 
0,1,... ,m with p m (q) = 1 and m > 0. 

Denote by Sk(q), k — 1,... ,m, the roots of p(q, •) and by a : Q —* M (or a p if it is necessary 
to clarify the dependence on the polynomial) the abscissa map of p, i.e. the maximal real 
part of the roots: 

a(q) := max 5R(s fc (g)), q G Q. 

k=l,...,m 

Equivalently, with i = ^/—T and s = x + iy write 
p(q, s) = Px{q, x , y) + ipct(q, x, y) 

for two real polynomials p^,p<^ € M[g,x,y] of total degree m. Then 

a(q) = max {a; 61:3 y eR: p^{q, x, y) = p^(q, x, y) = 0}, q € Q. 

We observe that function a : Q —» M is semi-algebraic and we define the basic closed 
semi-algebraic set 

z := {(g, X, y) e R n x M 2 : q G Q, p x (q, x, y) = p%(q, x, y) = 0}. 
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Remark 1. Set Z is compact, since Q is compact and p is monic in s. 

Now we can write the abscissa map as 

a(q) = max{x G M : 3y G M : (q,x,y) G Z}, q G Q. 

Since p is monic, its abscissa a is continuous, though in general not Lipschitz continuous. 
For example, for n — 1 and p(q, s) = s 6 + q the map a(q) is only Holder continuous with 
exponent | for small q. To be precise, a is always Holder continuous by the Lojasiewicz 
inequality [3], since Q is compact. 


3 Upper abscissa approximation 

3.1 Primal and dual formulation 

Given a polynomial p defined as above, the solution to the following linear programming 
(LP) problem gives an upper approximation of its abscissa function a on Q: 

p — inf / v(q)dq (1) 

veV(Q)jQ 

s.t. v(q) — x > 0 for all (g, x,y) G Z 

with ^(Q) denoting the space of continuous functions from Q to M. 

Remark 2. Since the continuous functions defined on compact set Q can be approximated 
uniformly by polynomials by the Stone-Weierstrafi theorem [II] §16.4.3], we can replace 
^(Q) in problem ([!]) by the ring of polynomials M[g]. 

The LP dual to problem (JTJ) reads 

P* = sup / x dp(q,x,y) (2) 

s.t. f q a dp = I q a dq , for all a G LP, 

Jz Jq 

where q a stands for the monomial q^q^ 2 • • • g“ n and ^ + (Z) denotes the cone of non¬ 
negative Borel measures supported on Z , identified with the set of all non-negative 
continuous linear functionals acting on ^ + (Z), the cone of non-negative continuous 
functions supported on Z. 

Remark 3. The constraint J z q a dp = Jq q a dq for all a G N n implies that the marginal of 
/i on Q is the Lebesgue measure on Q, i.e. for every g G ^(Q) it holds that 

/ 9(q)dp(q,x,y ) = / g(q)dq. 

Jz Jq 

In particular this implies that ||/i|| = vol Q where vol(-) denotes the volume or Lebesgue 
measure. 
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Lemma 1. The supremum in LP ([2]) is attained, and there is no duality gap between LP 
|7]j and LP |Ij), i.e. p = p*. 


Proof. The set of feasible solutions for the dual LP (J2]) is a bounded subset of + (Z ) 
with Z compact and therefore it is weak-star compact. Since the objective function is 
linear, its supremum on this weak-star compact set is attained. For elementary background 
on weak-star topology, see e.g. [21 Chapter IV], 

To prove that there is no duality gap, we apply [;2j Theorem IV.7.2], For this purpose 
we introduce the notation used in [2] in this context. There, the primal and the dual are 
written in the following canonical form: 

P* = sup (x, c)i p= inf (b, y) 2 

xeEi ye^2 

s.t. Ax = b, x e Ef s.t. A*y — c 6 Ff 

So we set Ei := ^(Z) with its cone Ef := ^ + {Z ). Then their (pre-)duals are Fi := c <o(Z) 
and Ff := ^ + (Z ) respectively. Similarly, we define E 2 := ~#(<2) and F 2 := ^(Q). 

Setting x := ^ G Fi, c := i 6 Fi, b 6 P 2 the Lebesgue measure on Q and y := v G F 2 , 
the linear operator A : Ei — > E 2 is given by x t—)■ 7Tgx where ttq denotes the projection 
onto Q, i.e., Ax(R) = x(£> x M 2 ) for all B e £>(Q). 

According to [21 Theorem IV.7.2] the duality gap is zero if the cone {(Ax, (x, c)i) : x 6 Ef } 
is closed in E 2 x M. This holds in our setup since x i —> Ax and x t—>■ (x, c)i are continuous 
linear maps and Ef = ^ + {Z) is weak-star closed due to the compactness of Z. So if 
for some a G E 2 , Ax n —^ a as n —> oo then from the dehnition of the mapping A and 
as (x n ) C Ef, one has ||x n || — > ||a|| as n —>■ oo (see Remark [ 3 ]). Therefore the sequence 
(x n ) C Ef is bounded and by Banach-Alaoglu’s theorem BE], it contains a subsequence 
(x llk ) C Ef that converges to some x e Ef for the weak-star topology. By continuity of 
the mappings A and c, the result follows. □ 

Remark 4. The infimum in LP (0 is not necessarily attained, since the set of feasible 
solutions is not compact. It is neither attained when we replace 1f(Q) with M[g], since a 
is non-Lipschitz, so in particular not polynomial. 

However, the infimum is attained if we replace 1f(Q) with M[g]d for d finite. Then, with 
M := miriggQ a(q ) > — 00 and v(q ) := v(q) — M we can rewrite LP ([Tj) as the equivalent 
problem 

inf / v(q)dq s.t. v(q) + M — x > 0 on Z. 

eeM [ q ] d jQ 

Now, any feasible v is non-negative on Q, so fQv(q)dq = ||€i|| Z/ i > 0 and for every Uel 
the set {u G M[g]d : R. > Jq v(q)dq and v(q) + M — x > 0 on Z } is closed and bounded in 
the strong topology, thus compact. Besides, due to the continuity of a, there always exists 
an R < 00 such that the mentioned set is not empty, hence the infimum is attained. 
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3.2 SDP hierarchy 


Let- do € N be sufficiently large. As presented in [9j, we can write a hierarchy of finite¬ 
dimensional convex semidefinite programming (SDP) problems for LP ([Tj) indexed by the 
parameter dsN, d > do'- 

Pd = inf / v d (q) dq 

V d ,CTo,(Tj,TsR,Tc5 Jq 

n 

s.t. v d (q) - x = a 0 (q,x,y) + a j(<L x ,y)( 1 - g 2 ) (3) 

3 =1 

+ Tn(q, X , y)ptR(q, x, y) + Tc>(g, x, y)p Q (q, x, y) 

for all (g, x, y) G M n x I 2 and with v d G R[g] 2 d, cr 0 e Tj[q,x,y\ 2d , e £[g, x, y] 2 d -2 for 
i = 1, • ■ ■, n and r R , r 0 G R[q, x, y] 2d - m . 

Remark 5. The quadratic module generated by the polynomials 1—gf,..., 1 — g 2 , ±psR, 
is archimedean by da Lemma 3.17], since it contains the polynomial f(q, x, y) := J2j= i(l — 
qj) — and the set {(g, x, y) G M n x M 2 : /(g, x, y) > 0} is compact. By [9] Theorem 

4.1] , this implies that the hierarchy converges, i.e. limbec p d = p- 

Remark 6. Note that SDP (|3]) is not equivalent to LP ([TJ) , not even with *if(Q) replaced 
by M[g] or M[g] 2 d in the latter, but it is a strengthening of it, meaning p d > p. To be more 
specific, SDP (J3]) is a reinforcement of the following LP: 

inf / v{q) dq s.t. v(q) -x>0 for all (g, x, y) G Z. 

^6K[<?]2d JQ. 

where we exchanged non-negativity for a specific certificate of positivity. See [9J Chapter 

4.2] for details. 

Example 1. The inffinum in SDP (|3]) is not necessarily attained, e.g. consider the 
polynomial p(q, s) =s 2 . Then p^(g, x, y) = x 2 — y 2 , p$s{q,x,y) — 2xy and Z = Qx{(0, 0)}. 
Obviously, the optimal solution to LP ([I| is v = 0. For SDP (j3j) we would want 

v(q) -x = cr 0 (g, x, y) + <7i(g, x, y)(l - g 2 ) + r^{q, x, y)(x 2 - y 2 ) + 2r a (g, x, y)xy, 

meaning {) = v = x + ao + (Ti(l — k 2 ) + t^x 2 — r^y 2 + 2r^xy with (To, sums of squares. 
This is impossible, since it would require the construction of the term — x which in this 
case is only possible as a summand of cr 0 . Then however we would always also produce a 
constant positive term. Practically this means that the multipliers cr 0 , oq, tsr, Tq blow up. 


Hence, an optimal solution might not exist, but we always have a near optimal solution. 
This means we should allow solutions v d with J Q v d (q) dq < p d + e.g. in the above 
example we would search for v = e for an e > 0 sufficiently small. 

Remark 7. The existence of an optimal solution depends on further conditions, like the 
ideal generated by the polynomials 1 — q 2 , p^ and pq being radical, and goes beyond the 
scope of this paper. The interested reader is referred to the proof of [7, Lemma 1] for 
further details. 
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Ill the following theorem we prove that the associated sequence of solutions converges: 

Theorem 1. Let Vd G R[g]2d be a near optimal solution for SDP (J3]) , i.e. f Q Vd(q) dq < 
Pd+ and consider the associated sequence ( Vd)d>d 0 C L 1 (Q). Then tq converges to the 
abscissa a in L 1 norm on Q. 


Proof. Recall that p* = p according to Lemma [lj First we show that p — f Q a(q) dq. For 
every (q, x,y) G Z we have x < a(q) and since f z q a dp = f Q q a dq for all a G N n which 
means that the marginal of p on Q is the Lebesgue measure on Q (see Remark [3j), it 
follows that for every feasible solution p G A% + (Z) it holds that 

x dp(q, x, y) < / a(q) dp(q,x,y) = / a(q)dq. 

Jz Jq 

Hence p < f Q a(q) dq. On the other hand, for every q G Q there exists {q,x q ,y q ) G Z such 
that a(q) = x q . Let p* be the Borel measure concentrated on (q,x q ,y q ) for all q G Q, i.e. 
for A in the Borel sigma algebra of Z it holds 



p*(A) := l A (q,x q ,y q ). 

Then p* is feasible for problem (J2| with value 

/ x dp*(q, x, y) = / a(q) dq, 

Jz Jq 

which proves that p > f Q a(q ) dq, hence p — f Q a(q) dq. 


Next we show convergence in L l . Since the abscissa a is continuous on the compact set Q, 
by the Stone-Weierstrafi theorem [TTJ, §16.4.3] it holds that for every e > 0 there exists a 
polynomial h £ G R [q] such that 

sup \h e (q) ~a(q) \ < 

q&Q ^ 

Hence, the polynomial v £ := h £ +e satisfies v e — a. > 0 on Q and we have v e (q) — x > 0 on 
Z. Since the corresponding quadratic module is archimedean (see Remark [5]), by Putinar’s 
Positivstellensatz P Theorem 2.5] there exist <7q, cb G Z[q,x,y\, r^, r§ G R[g, x,y] such 
that for all ( q , x, y ) G R n x R 2 we can write 

n 

v £ (q) - x = a £ 0 (q, x, y) + ^ ^ e (g, x, y)( 1 - q]) 

j =i 

+ x, y)ptR(q, x, y) + r|(g, x, y)p$(q, x, y). 

Therefore, for d > d £ [ dc ^ t,£ ] the tuple (v e , erg, cr*-, r^, r|) is a feasible solution for SDP 
([3]) satisfying 

0 < [ (v £ (q) - a(q)) dq<^ f dq. 

Jq 2 Jq 
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Together with f Q a(q) dq = p < p d which is due to the first part of the proof and p d being 
a strengthening of p, it follows that whenever d > d £ it holds that 

0 < pd~ [ a(q ) dq< [ (v e (q) - a{q)) dq<^ [ dq. 

JQ Jq 2 Jq 

As £ > 0 was arbitrary, we obtain lim^oo p d = Jq a(q ) dq and since a <v d for all d, this is 
the same as convergence in L 1 : 


0 < lim ||nd — a\ 

d—>oc 


i = lim 

d—>oo . 


'Q 


I v d (q) ~ a (q)\ dq 


= lim / (v d (q) 

a—>-oo JQ 


a(q)) dq < lim (p d + 

d->cxD v (i 


IQ 


a(q) dq = 0. 


□ 


For linear systems, a polynomial is called stable if all its roots lie in the open left part of the 
complex plane, i.e. if its abscissa is negative. Hence for a polynomial with parameterized 
coefficients, as we consider in this paper, the stability region is the set of parameters for 
which the abscissa is negative, in our notation 

{q E Q: a(q) < 0}. 

The following statement on polynomial inner approximations of the zero sublevel set of 
the abscissa function follows immediately from the L 1 convergence result of Theorem [lj 
see also [7]. 

Corollary 1. Let v d e M[g] 2(i denote, as in Theorem [7], a near optimal solution for SDP 
©■ Then {q G Q : v d (q) < 0} C {q G Q : a{q) < 0} and lim^^ vol {q <G Q : v d (q) < 
0} = vol {q G Q : a(q) < 0}. 


3.3 Examples 

As stated in Corollary [lj while approximating the abscissa function from above we also 
get an inner approximation of the stability region. The authors of |7] surveyed a different 
approach. They described the stability region via the eigenvalues of the Hermite matrix of 
the polynomial and approximated it using an SDP hierarchy. In the following examples we 
compare the two different methods and highlight the specific advantages of our abscissa 
approximation. 

In this section and in the remainder of the paper, all examples are modelled by Yalmip 
and solved by Mosek 7 under the Matlab environment, unless indicated otherwise. 

Example 2 (The damped oscillator [6]). Consider the second degree polynomial depending 
on n — 1 parameter q G Q — [—1,1]: 

p : s i —y p(q, s) = s 2 + 2 qs + 1 — 2 q. 
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Figure 1: Abscissa (black) and its polynomial upper approximations of degree 4 (left, gray) 
and 10 (right, gray) for Example [2] The quality of the approximation deterioriates near 
the minimum, where the abscissa is not Lipschitz. 


Then Z = {(g, x, y) G [—1,1] x M 2 : x 2 — y 2 + 2 qx + 1 — 2 q = 2 xy + 2 qy = 0} and the 
corresponding hierarchy of SDP problems (J3]) reads 

Pd = inf f v d (q) dq 

Vd,VO,<Tl,TSR,TQ J_l 

s.t. v d (q) - x = cr 0 (q,x,y ) + cq(g, x, y)(l - q 2 ) 

+ m(q, x, y)(x 2 - y 2 + 2qx + 1 - 2 q) + r 0 (g, x, y)(2xy + 2 qy) 


for all ( q,x,y) G M 3 and with v d G M[g] 2 d, cro e T,[q,x,y\ 2 d, cri 6 E[g,x,y] 2 rf -2 and 
' r 5 R>' r o ^ M[g, x, 2/] 2 rf— 2 - Apart from that, we only need the moments of the Lebesgue 
measure on [—1,1] for a successful implementation. These are readily given by 





l - (-i )" +1 

ot + 1 


meaning that f\ v d (q) dq = J2a =l v da z a with v da denoting the coefficient of the monomial 
q a of v d . See Figure[l]for the graphs of the degree 4 and 10 polynomial upper approximations 
of the abscissa. 


For the Hermite approximation we compute the Hermite matrix H of p (see fH] for details) 
H{q > ~ { 0 4 q 

and write the hierarchy of optimization problems as presented in [7]: 
max f g d (q) dq 

9d& 0i cr l,' r J — 1 

s.t. u T H(q)u - g d (q) = a 0 (q, u ) + a^q, n)(l - q 2 ) + r(q, u)( 1 - u T u) 
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for all (q,u) £ [-1,1] x M 2 and with g d £ M[g] 2 d, cr 0 £ 'E[q,u\ 2 d, oq e T 1 [q,u } 2d -2 and 
r G R[q,u] 2 d- 2 - Already for d = 6 we observe a close match between the genuine stability 
region, which is {q £ [—1,1] : a(q) < 0} = (0,1), the Hermite inner approximation {q £ 
[—1,1] : —ge(q) < 0}, and the polynomial upper approximation {q £ [—1,1] : Vio(d) < 0}. 
These three intervals are visually indistinguishable, so we do not represent them graphically. 




Figure 2: Abscissa (black) and its polynomial upper approximations of degree 6 (gray, 
left) and 12 (gray, right) for Example [3] The quality of the approximation deterioriates 
near the points of non-differentiability of the abscissa. 

Example 3. Consider the polynomial 

p : s (->• p(q, s) = s 3 + 1 s 2 + q 2 s + (q - \)q(q + §) 

for q £ Q = [—1,1], The abscissa function a(q) of p is not differentiable at three points 
and therefore it is rather hard to approximate in their neighborhoods. In Figure [2] we see 
the abscissa and its polynomial upper approximations of degrees 6 and 12. Comparing 
the genuine stability region {q £ [—1,1] : a(q) < 0}, the polynomial inner approximation 
{q £ [—1,1] : v 12 {q) < 0} and the Hermite inner approximation {q £ [—1,1] : — gio(q) < 0}, 
we observe, maybe surprisingly, that the approximations are very similar and miss the 
same parts of the stability region. These are not reproduced graphically. 

Remark 8. Evidently, the approach via the Hermite matrix does not tell us anything 
about the abscissa function itself besides from where it is negative. As an illustration 
consider a polynomial of the form p(q, s) = s 2 + po(q) for n — 1. Then p(q, •) has either 
0 as a multiple root, two real roots (of which one is positive) or only imaginary roots, 
i.e. the stability region of p is empty and its Hermite matrix H(q) is zero. Therefore the 
eigenvalues and their approximation g d are also zero for every d. In contrast, the upper 
abscissa approximation v d gives a suitable approximation for the abscissa function. 

On the other hand, practical experiments (not reported here) reveal that computing 
the abscissa approximation is typically more challenging numerically than computing 
the Hermite approximation. For instance, computing the upper abscissa approximation 
may fail for polynomials with large coefficients, while the Hermite approximation keeps 
providing a proper inner approximation of the stability region. 
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Figure 3: Abscissa (dark, below) and its polynomial upper approximations of degrees 
6 (left, transparent) and 10 (right, transparent) for Example [4] We observe that the 
approximation deteriorates near the regions of non-differentiability of the abscissa. 


Example 4. Consider the polynomial 

p: p(q , s) = s 3 + {qi + §)s 2 + q\s + q x q 2 

depending on n = 2 parameters q G Q — [—1, l] 2 . Then Z = {(q,x,y) G [— 1,1] 2 xl 2 : 
x 3 - 3 xy 2 + (qi + |)a; 2 - (q\ + § )y 2 + qfx + gig 2 = -y 3 + 3 x 2 y + 2(gi + | )xy + qjy = 0}. 
In Figure [3] we represent the graphs of the abscissa a and its polynomial approximations 
Vq and r? 10 . In Figure [4] we represent the stabilizability region, i.e. the zero sublevel set of 
the abscissa {q G [—1, l] 2 : a(q) < 0} (dark gray region), the degree 8 Hermite sublevel set 
{q G [—1, l] 2 : —gs(q) < 0} (light gray region, left) and the degree 10 polynomial sublevel 
set {q G [—1, l] 2 : Vio(q) < 0} (light gray region, right). 

Remark 9. In the examples we always chose lower degrees for the Hermite approximation 
than for the upper abscissa approximation. The Hermite approximation converges relatively 
fast making it unnecessary to consider higher degrees, especially since they require much 
more time. On the contrary, the upper abscissa approximation usually needs higher degrees 
to provide a useful approximation, but it is faster to compute. 
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Figure 4: Stabilizability region (dark gray region) and its inner approximations with degree 
8 Hermite (light gray region, left) and degree 10 upper polynomial approximation (light 
gray region, right). 

4 Lower abscissa approximation 


At first thought, finding a lower approximation for the abscissa map might sound like a 
straightforward task, since one is tempted to just solve the analogue of LP (JTj): 

sup / w(q)dq (4) 

weV(Q) JQ 

s.t. x — w(q ) > 0 for all (g, x, y) G Z. 


This, indeed, gives a valid lower bound on the abscissa function, however in general a very 
bad one since it is not approximating the abscissa but the minimal real part of the roots 
of p. To understand the reason we recall that 

Z = {(<?, x, y) G M n X M 2 : q G Q, p x (q, x, y) = p^(q, x, y) = 0} 

and therefore this set contains all roots of p and not only those with maximal real part. 

Example 5. On the left of Figure [5] we show the degree 12 solution to the SDP hierarchy 
corresponding to LP Q for the polynomial p(q, s) = s 2 + 2 qs + 1 — 2q of Example [2j 
which gives a tight lower approximation to the abscissa only in the left part of the domain, 
corresponding to a pair of complex conjugate roots. We observe that the SDP solver 
Mosek does not return a correct answer for this particular problem, and we had to use 
the SDP solver SeDuMi instead in this case. On the right of Figure [5] we show the 
degree 12 solution to the SDP hierarchy corresponding to LP Q for the polynomial 
p(q, s) = s 3 + |s 2 + q 2 s + (q — \)q{q + of Example [3j The lower approximation is 
nowhere tight, due to the presence of roots with real parts smaller than the abscissa. 
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Figure 5: Real parts of the roots (black) and degree 12 polynomial lower approximations 
(gray) for the second degree polynomial (left) and third degree polynomial (right) of 
Example [5} 


To find a tighter approximation for the abscissa map from below we pursue two different 
approaches: 


• First, we reformulate the set Z with the help of elementary symmetric functions, in 
order to have access to the roots directly. This is a very neat way with options for 
variation, such as approximating the second largest real part of the roots from above 
or below, but it also includes many additional variables and it is therefore not very 
efficient when implemented. However, it can be useful for small problems. 

• Second, we restrict LP Q further using the GauB-Lucas theorem, i.e. instead of 
Z we use a subset of Z which contains only the roots with the abscissa as its real 
parts. This approach is much more complicated, relies on assumptions and one 
needs to solve two optimization problems in order to get the lower approximation. 
Nevertheless, the implementation is much faster and it can also be used for bigger 
problems. 


4.1 Lower approximation via elementary symmetric functions 

4.1.1 Problem formulation 


Let us derive another description of the set of roots of p which allows us to pick single 
roots according to the size of their real part. For this purpose let us recall the definition 
of our polynomial: 

m 

P ■ s p(q, s ) := J2 Pk{<l)s k with p m (q) = 1. 

k =0 
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Following the notation of the previous chapters, we denote the roots of p(q, •) by s k {q), k = 
1 ,,m and split them up into their real and imaginary parts, s k (q) = x k (q) + iy k (q) with 
x k(q),Vk(q) £ ®h To simplify notations we omit the dependence on q whenever it is clear 
and write only s k , x k and y k . 

Now we write the coefficients of the polynomial as elementary symmetric functions of its 
roots: 


Pm-k(q) = (-l) fc s h s i2--' s ik’ k = l,...m. 

This allows us to define the set of roots of p in the following way, where we can order the 
roots arbitrarily: 

Z'o := {(q, X!, ..., x m , yi, ..., y m ) e Q x R m x R m : x k < x m , k = 1,..., m - 1, 

Pm-k(q) = (~ l) fc X! s h s h-■ ■ s i k i k = l,...m}. 

To avoid complex variables si k in the description of the set, we replace them by si k = 
xi k + iyi k and split up the sum Xu <h<-<i k <m s h s h ''' s h i n rea ^ an( i imaginary parts. 
The latter would be zero, since all p m _fc(g) are real. In the sequel we omit this procedure, 
since it would only complicate the notations. 

For illustrative reasons let us fix q for a moment. Then the set Z' 0 contains only one 
element (q, aq,..., x m , yi ,..., y m ). For this it holds that x rn = a(q) and the points 
(q, x k , Vk)-, k = 1,... , m, are exactly the elements of Z. 

Remark 10. One could order the roots further by adding more conditions, like for example 
Xk < x m ~ i, k = 1, • • • pm —2. Then one could also access the root with the second largest 
real part. Of course, this would imply another rn — 2 constraints in an implementation 
and therefore this would slow down further the solving process. 


In theory, m variables suffice to characterize the roots of a real polynomial via the 


elementary symmetric functions, but since we need a 
identify the maximal one, we can only eliminate 
We set 


1 variables x k explicitly in order to 
:= maxjr e Z | r < variables. 


Uk-i 



2.. .. pm 

2.. .., m - 3 and y m _ 2 


if m is even 
-ym -i - Vm if rn is odd, 


meaning we decide which roots will be pairs in case they are complex. Note that it is 
necessary to keep y m , since we defined x m as the abscissa and we do not know whether s m 
is real or not. In fact s m (q) can be real for some q and complex for others. 


Remark 11. 
eliminate 


Iven though we know for m odd that one root must be real, we cannot 
variables, since it might happen that s m (q) is the single real root for some 


q while it is complex for other q. 
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Now we can write the set of roots with less variables and less constraints. As above we 
keep the variables Sk in the description of the set for readability reasons, but remark that 
with the reduced amount of y variables the constraints 0 = ^s(J2i<i 1< ... < i k < rn Sj 1 Sj 2 • • • Si k ) 


for k — 1,..., yj are superfluous. We have 


rm-i 

Z 0 :^{(q,x 1 ,.. . ,x m ,y 2 ,V 4 , ■ ■ ■ ,y 2 ^m^y m ) G Q x R m x Ml 2 I ; 

Xk < x m , k = 1 ,... ,m - 1 , 

Pm-k(q) = (~ l) fc 51 S h Si 2 ■ ■ ■ Si h , k = l,... m}. 


Example 6. For m — 3 the set Z a is given by 


Z a ={(q,x 1 ,x 2 ,x 3 ,y 2 ,y 3 ) £ Q x R 3 x i 2 : ii < i 3 , x 2 < x 3 , 

~P 2 (q) = xi + x 2 + x 3 , 

Pi(q) = xix 2 + x 3 x 3 + x 2 x 3 + y 2 + 2 / 22/3 + y \, 

- p 0 (q) = xxx 2 x 3 + {-Xi + x 2 + x 3 )y 2 y 3 + x 2 y\ + x 3 y\, 
0 = (xi - x 2 )y 2 + (xi - x 3 )y 3 , 

0 = (a?i — x 2 )x 3 y 2 + (xi - x 3 )x 2 y 3 + y\y 3 + y 2 y\}. 


To clarify the formula also for m even, we write Z 0 down explicitly for m = 4: 

={{q,x 1 ,x 2 ,x 3 ,x 4 ,y 2 ,y 4 ) G Q x M 4 x M 2 : xi < x 4 , x 2 < x 4 , x 3 < x 4 , 

~ P:i(q) = X! + x 2 + x 3 + x 4 , 

p 2 (q ) = £ 10:2 + £ 1^3 + X 1 X 4 + x 2 x 3 + 2:2X4 + X3X4 + y\ + 7/4, 

- pi(g) = xix 2 (x 3 + x 4 ) + (xi + x 2 )x 3 x 4 + (xi + x 2 )y 2 A + (x 3 + x 4 )y 2 , 
Po(q) = X!X 2 X 3 X 4 + (xi - x 2 )(x 4 - X 3 )y 2 y 4 + x 4 x 2 y 2 + x 3 x 4 y 2 + y 2 y 4 , 
0 = (xi - x 2 )(x 3 + x 4 )y 2 + (xi + x 2 )(x 3 - x 4 )y 4 , 

0 = (xi - x 2 )(x 3 x 4 + 2/4)2/2 + ( x 3 - x 4 )(xix 2 + 2/2)2/4}- 


Here we have set 2/1 = — 2/2 and 2/3 = — 2 / 4 , so the constraint 0 = 3f(Xa<j 1 <...<j fc < m s/ 1 sj 2 • • • sj fc ) 
for k — 1 is obviously superfluous, because it reduces to 0 = 0. The second superfluous 
constraint is the one for k = 2, that is 0 = (x 4 — x 2 )y 2 + (x 3 — x 4 )y 4 , since we have x\ = x 2 , 
respectively x 3 = x 4 , in the case s 2 , respectively s 4 , is complex. 


Finally, we can reformulate LP Q in such a way that it provides a proper approximation 
of the abscissa function from below: 

p = sup / w(q) dq (5) 

w&V(Q)-JQ 

s.t. x m - w(q) > 0 for all (q,x 1 ,..., x m , y 2 , y 4 ,.. .,y 2 yrn_^y m ) e Z a . 

With the notation of Section l3Tl its dual LP reads 

P bif / x m dp(q, Xi ,..., x m , 2 / 25 1/4; ■ • •; y 2 \ m 1 ; Pm) (6) 

fi€.//+(Zo) Jz o L 2 J 

s.t. f q a dp = f q a dq, for all a G N n . 

Jz 0 Jq 
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In analogy with the upper approximation we have no duality gap and the infimum is 
attained: 

Lemma 2. The infimum in LP is attained, and there is no duality gap between LP 
and LP |6j) ; i.e. p = p*. 

Since Z 0 is compact, the proof is identical to that of Lemma [lj 

Remark 12. For the same reasons as for the upper approximation ([Tj) , the supremum in 

8 is not attained for or M[g], but it is attained for M.[q\ d with d finite. See Remark 

vith M := min, /e Q a(q) — N for an IV 6 N sufficiently large, and R := f Q ( a(q ) — M) dq. 


4.1.2 SDP hierarchy 


Let do G N be sufficiently large. Then for d G N, d > d 0 the corresponding hierarchy of 
SDP problems reads 


Pd = sup / w d (q) dq 

j °7 c i&xi ^ Q. 


(7) 


m— 1 


S.t. X m - w d (q ) = cr 0 + °'j( 1 - Tj) + E °x k {Xm ~ X k ) 


+ E (-1 ) k pm-k{q) ~ 3? I E 


s li s h ' ' ' s lk 


k =1 


E 

fc= LfJ 


E 

i l<li<l2<--<l k <m 


l<h<h<-<lk<m 
Sli^h ' ' ' &l k 


rrn i 

for all (q,Xi,...,Xm,y 2 ,yi, ■ ■ ■ ,2/ 2 [^J j 2/m) G 1" x I m x Ml 2 and with w d G it^[yj 2 d, 

0-0, Trfc ^ ^[*?) ^T, ■ ■ ■ ■ X m , 2/2, 2/4) • • • ) 2/m]2d for k 1, . . . , 7T7. 1, CTj. G . . . , X ?n , 

2/2,2/4) • • • ,ymhd -2 for k = r R , k G K[g, aq,..., x m , y 2 , 2/4, • ■ ■, Vmhd-k for k = 

1,..., m and for k = yj ,..., m. 

Remark 13. As in Remark [6j SDP (|T]) is a strengthening of LP ([5]), meaning p d < p. 
Also as in Remark [5| the quadratic module corresponding to Z Q is archimedean, i.e. 
ffi 11/ >oc Pd P- 


We conclude the section with the following result: 

Theorem 2. Let w d G M[g] 2 d be a near optimal solution for SDP ([7]) ; i.e. fQW d (q) dq > 
Pd — h and consider the associated sequence (w d ) d > do C L 1 (Q). Then w d converges to a in 
L 1 norm in Q. 


Unsurprisingly, one can prove this result in exactly the same way as Theorem [lj so we do 
not detail the proof here. Remark that the first part of the proof can be shortened, since 
fz 0 dp(q, Xi, ..., x m , ?/2, 2/4, • • •, 2/m) f z u ^Pi.Q .1 Xi, ■ ■ ■, x m , |/ 2 , 2 / 4 , • • • , 2/m)- 
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4.1.3 Examples 


Just as the upper abscissa approximation automatically approximates the stability region 
from inside, the lower approximation gives, as a side effect, an outer approximation. In 
this section we will examine similar examples as for the upper approximation. 



q 



Figure 6: Abscissa (black) and its polynomial lower approximations of degree 6 (gray, left) 
and 10 (gray, right) for Example [7j The quality of the approximation deteriorates near 
the minimum, where the abscissa is not Lipschitz, compare with Figure [l] 

Example 7. As in Example [2] consider the polynomial 
p : s i —^ p(q, s) = s 2 T 2 qs + 1 — 2 q. 

We have t/i = -t/ 2 , so := {(g, x \, x 2 , y 2 ) G Q x M 3 : x\ < x 2 , -2 q = x± + x 2 , 1-2 q — 
xix 2 +y 2 i 0 = (xj — x 2 )y 2 }. In Figure [6] we see the graphs of the degree 6 and 10 polynomial 
lower approximations obtained by solving SDP ([7]). As in Example [5| we observe that the 
SDP solver Mosek does not return a correct degree 10 polynomial, and we had to use the 
SDP solver SeDuMi instead in this case. Due to the rather big amount of variables and 
constraints, computing the degree 10 solution is already relatively expensive, with a few 
seconds of CPU time. 

Example 8. As in Example [3] consider the polynomial 

p : s h+ p(q, s) = s 3 + |s 2 + q 2 s + (q — \)q{q + |). 

With yi = — 7/2 — 1/3 we calculate Z Q as in Example [6} In Figure [7] we see the graphs of 
the degree 6 and 10 polynomial lower approximations obtained by solving SDP (JtJ) . The 
computation time to get the degree 10 solution is around 15 minutes, which is arguably 
not a good compromise given the quality of the approximation. 

Remark 14. As for the upper abscissa approximation, we observe practically that the 
implementation for the lower approximation is rather sensitive to polynomials with large 
coefficients. 
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Figure 7: Abscissa (black) and its polynomial lower approximations of degree 6 (gray, left) 
and 10 (gray, right) for Example [8j The quality of the approximation deteriorates near 
the minimum, where the abscissa is not differentiable, compare with Figure [2j 

Example 9. As in Example [4| consider the polynomial 

P ■ S'-i p{q, s) = s 3 + (q 1 + §)s 2 + g 2 s + q\q 2 . 

Since we have m — 3, the set Z Q is again given in Example [6} In Figure [8] we see the outer 
approximation of degrees 6 and 10 obtained by solving SDP (j7j). We notice an opening 
in the approximation of the stability region in the lower half of the picture. This is due 
to a being zero and non-smooth for q\ = 0, meaning a(0, qf) = 0. This phenomenon also 
incapacitates w$ to get tighter to a for q 2 > 0 than we observe in the upper half of the 
picture. 


4.2 Lower approximation via Gaufi-Lucas 

4.2.1 Problem formulation 


As indicated above, we want to find a semi-algebraic subset of Z which contains only 
those roots of p whose real part is maximal. This means that, in contrast to the approach 
of Section 4.1, we will not rephrase Z, but formulate further constraints. 


In order to do this we must distinguish between the roots of p(q, •) according to the size of 
their real parts. For this purpose we use the following result: 

Theorem 3 (GauB-Lucas). The critical points of a non-constant polynomial lie in the 
convex hull of its roots. 


We refer to [5] for further information and a proof. Let us denote the derivative of p(q, s) 
with respect to s by p'(q,s). By Theorem |3j the roots of p'(q,-) are contained in the 
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Figure 8 : Stabilizability region (dark gray region) and its degree 6 outer approximation 
(light gray region, left) and degree 10 outer approximation (light gray region, right). 
Compare with Figure [4} 

convex hull of the roots of p(q, •). It follows readily that the abscissa cy of p' lies below 
the abscissa a p of p: 

a p >(q) < a p (q) for all q G Q. 

However, p may have some roots with real part strictly smaller than a p and strictly bigger 
than dpi, meaning that the root whose real part is the abscissa is not the only one whose 
real part lies above cy. Of course, this cannot happen for real polynomials R-)-R because 
of monotonicity, and neither for complex polynomials of degree 2. But, for example, for 
n = 1 the polynomial p(q, s) = s 4 + ( q 2 + l)s + q has two roots with different real parts 
greater than ay for q G [—1, —0.4], 

To prevent the lower abscissa approximation from converging to the real part of a root 
smaller than the abscissa, we make the following assumption: 

Assumption 1. None of the real parts of any root of p lies strictly between a p and cy, 
i.e. x }a p fq),a p (q)[ for all ( q,x,y ) G Z. 

Remark 15. Unfortunately, we do not know how restrictive this assumption is. For n = 1 
it was rather difficult to find examples that violate it. 

Now let v G ^(Q) be a near optimal solution to LP (|TJ) for the polynomial p\ meaning 
f Q v(q) dq < p + e for an e > 0. Then, v is an upper approximation of the abscissa ay of 
pi. We define the following subset: 

Z r := {(< 7 , x, y) G Z \ x - v{q) > 0}. 

In order to see where we are going, let us pretend for a moment that we have an optimal 
solution. Then, under Assumption [lj the set Z r would contain exactly the points (q, a(q),y q ) 
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with y q denoting the imaginary part of the root of p(q, ■) with maximal real part. Hence, 
the solution to the following LP would give a lower approximation of the abscissa function 
a p of p: 

sup / w(q) dq (8) 

weV[g] J 2 

s.t. x — w(q ) > 0 for all (q, x, y ) G Z r . 

Since v might not be optimal, the projection of Z r onto Q can have holes of volume e. As 
a consequence, w might not be a valid lower bound of the abscissa on these holes. 

Taking this into account, we build an SDP hierarchy for LP (|8]) in the next section. The 
issue is that we have to consider the hierarchy for the upper approximation of Gy first and 
the solution to it might interfere with a p . 


4.2.2 SDP hierarchy 


For d' 0 G N sufficiently large we denote by v d >, d’ > d’ 0 , the solutions to SDP (J3]) for the 
polynomial p'. Thus, the v d > are polynomials in and by Theorem [l] the sequence 

(vd')d’e N converges to Gy from above in L 1 norm. 

Next, we want to describe the set Z r via the polynomials v d ' in order to have an irnple- 
mentable problem, i.e. we define 

Z r , d ' '■= {(q,x,y) G Z : x - v d >(q) > 0}. 

Of course, the set Z rd ; is highly dependent on the quality of v d / and hence on the choice 
of d!. Evidently, Z r d i is a subset of Z ri possibly strictly. To ensure that Z r d / contains all 
roots of p with the abscissa as their real parts we need v d > < a p . However, in practice this 
is impossible in some cases: 

Example 10. The abscissa a p of p{q,s) = (s 3 + q) 2 and the abscissa ay of p' coincide 
and have a point of non-differentiability at q — 0. As another example consider the 
polynomial p(q, s) = s 4 + qs for which both a p and ay are not differentiable at q = 0 and 
ey(0) = cy(0) = 0. 


For these examples we cannot achieve v d ' A cl p with d' finite, since v d > is a polynomial and 
therefore differentiable everywhere. 

As a consequence, we formulate another assumption. In general, the points that may cause 
problems are the ones where a p and ay coincide, i.e. the points of the set 

V ■= {q G Q | a p (q) = a p ,(q)}. 

On this set the polynomial v d ' should approximate ay perfectly for a finite d'. Calling a 
solution v d t near optimal if it satisfies fgVd'iq) dq < p d > + Jy, we assume: 
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Assumption 2. There is a near optimal solution v d i to SDP ([3]) for the polynomial p' 
with d! finite such that v d i and a p > coincide on V. 

Remark 16. A sufficient condition for a violation of Assumption [2] is the existence of a 
value of q for which cy is not differentiable and a p (q) = a p fq ). This is the case for the 
examples given above. Note also that they are of degenerate nature. 

To face another issue, we denote the projection of Z r ^ onto the set Q by 7TQ(Z rid i), i.e. 

KQ{Zr,d') = {q e Q : 3x,y G M : (q,x,y) G Z r ^}. 

Since v d / converges to a in L 1 , but not necessarily uniformly, it might have spikes or 
similar irregularities, meaning that the set Q \ 7iQ(Z rid /) is not empty. However, the L 1 
convergence of v d i, or more precisely the convergence in measure, implies that there is a 
subsequence (v d ')ie n which converges to cy almost uniformly (see e.g. [TJ, Theorem 2.5.3]). 
In other words, for all 5 > 0, there exists a set As in the Borel sigma algebra of Q such 
that Sa s dq < 5 and v d > converges uniformly on A d to ay when l — » oo, where A d is the 
set-theoretic complement of As in Q. With this notation we have 

*Q(Zr,d ') C A° 5 C Q. 

Lemma 3. Let Assumption hold. Then , for every 5 > 0 there is a finite rf'eN such 
that v d r < a p on A d . 


Proof. Fix S > 0. Obviously we want 

0 < a p (q) - v d fq) = a p (q) - a p fq ) + a p fq) - v d >(q) (9) 

for every q G A^ C Q. By Theorem [3j we have a p (q) — a p >(q ) > 0 for all q G Q. Otherwise, 
the difference a p i (q) — v d > (q) is negative by construction, but due to Theorem [I] we find a 
subsequence v d > converging uniformly to a p i on Af. Hence, there is a finite d[, such that 
Q is fulfilled for all q G {q G A$ : a p {q ) > a p /(q)}. Because of Assumption [ 2 ] there is also 
a finite d' such that a p fq) — v d fq) vanishes on {q G Af : a p (q) = a p fq)}. Taking d' v > d' 
with l' > l* completes the proof. □ 

Remark 17. Choosing d' according to Lemma [ 3 ] implies ttq{Z t ^) = A^. 


Under Assumption [I] and [ 2 ] and with an appropriate choice of d! (depending on 5) the 
solution to the following LP gives a lower approximation of the abscissa function a p of p 
on the set A£ C Q: 


p d > = sup / w(q) dq (10) 

w£^[q] Jl TQ(A,d') 

s.t. x — w(q) > 0 for all ( q,x,y ) G Z T( y. 


Remark 18. Note that under Assumption [TJ LP (10) always provides a proper approximation 
for the abscissa a p from below on TiQ(Z rd /), but this might not be very useful, since for 
bad v d ' this set may have big holes or even be empty. To achieve suitable results on Afj 
we need Assumption [2] and an appropriate d', meaning a sufficiently good v d > ensuring 
KQ.{Z, rA ') = As- 
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Iii analogy with (|2]), the dual LP reads 

Pd'= J n i J xdp(q,x,y) (11) 

^JK+{Z r d ,) JZ r d i 

s.t. [ q a dp = I q a dq , for all a G N n 

•* Z T-,d' "' 7 r Q(- Z r,d') 

with the notation of Section 13.11 

Lemma 4. The infimum in LP 0 is attained, and there is no duality gap between LP 
(10) and LP 0 > «-e. p d , = p* d ,. 


Since Z rd i is a compact subset of Z , we can mimic the proof of Lemma [l] in order to get a 
proof of Lemma [4j 

Remark 19. As in Remark |4j the supremum in LP (10) is not attained for Q ) or 
M[g], but it is attained for M[g] rf with d finite. To adjust the proof of Remark |4j set 
M := miiiqgQ a(q) — N for an N G N sufficiently large, and R fa ( a(q) — M) dq as in 
Remark [12] 


Finally, for d! as in Lemma [3] and do > d' sufficiently large we can write an SDP hierarchy 
indexed by d G N, d > d^: 


Pd',d = sup / w d (q) dq 


( 12 ) 


s.t. x-w d (q) = a Q (q,x,y) + ^2<Tj(q, x ,y)(l - q]) + <Te(q,x,y)(x - v d '(q)) 

3 =i 

+ PR (q, x, y)p*t(q, x , y) + r Q (g, x, y)p^{q, x, y) 

for all ( q,x,y ) G M n x M 2 and with w d G M[g] 2d , cr 0 G S [q,x,y] 2d , crj G Z[q,x,y ] 2d -2 for 
j = 1,..., n, an G E[g, x, y] 2 d-d' and r», G M[g, x, y} 2 d-m- 


Remark 20. As in section 3.2 SDP (12) is a strengthening of LP (10), meaning p d i d < p d <. 
Besides, the archimedean quadratic module corresponding to the set Z is contained in 
the quadratic module corresponding to Z r ^. Hence, the latter is also archimedean, i.e. 
Bill,/ >oc p d ' d p d ' Pd 1 ' 

Remark 21. For numerical applications one can assume that As is empty and substitute 
Ka{Z rA ') by Q. 


The associated sequence converges: 

Theorem 4. Let Assumptions [7] hold and let Aff and d! be as in Lemma [5j Let 

w d G M[g] 2 d be a near optimal solution for SDP (12), i.e. f Q w d (q ) dq > p dd i — l ,. Consider 
the associated sequence ( w d ) d > do C L 1 (Q). Then w d is a valid lower bound of a p on Aff 
and it converges to a p in L 1 norm on Af. 


The proof of this result is very similar to the proof of Theorem [TJ so we omit it. Note that 
by Lemma [ 3 ] every feasible solution to SDP (12) is a valid lower bound of a p on Af and 
that we have 7iQ(Z rjd >) = Af due to our choice of d!. As for the proof of Theorem [ 2 J the 
first part can be shortened, since for every ( q,x,y ) G Z rd r it holds that x = a(q). 
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4.2.3 Examples 



q q 


Figure 9: Abscissa (black) and its polynomial lower approximations of degree 6 (gray, left) 
and 12 (gray, right) for Example 11 The quality of the approximation deteriorates near 


the minimum, where the abscissa is not Lipschitz, compare with Figures g and g. 
Example 11. As in Examples [2] and [7] consider 

p : p(q, s) = s 2 + 2 qs + 1 — 2 q. 


Assumption g is naturally fulfilled, since p is of degree 2. In the same way, Assumption 
gis fulfilled, since a p r(q) = -q is polynomial. We have Z r , d! {(+ *+ y) g [ 1 , l] x M . 
x — Vd'(q) >0, x 2 — y 2 + 2 qx +1 — 2 q — 2 xy + 2 qy = 0} and the corresponding SDP (g 
reads 

Pd',d = sup [ Wd(q) dq 

s.t. X - w d (q) = a 0 (q , x, y) + <+(+ x, y)( 1 - q 2 ) + a^(q, x, y){x - v d >{q)) 

+ r»(g, x , y)(x 2 - y 2 + 2 qx + 1 - 2 q) + r 0 (g, x, y)(2xy + 2 qy) 

for all ( q,x,y ) G M 3 and with G M[g] 2 d, G E[g,x, y\ 2d , G E[g,x,2/] 2 rf-2, Uc G 
E[g, x, 2/] 2 rf—rf 7 an d tsr, rcj G M[g, x, y] 2 d- 2 - Due to the simplicity of ey it suffices to choose 
d! = 2. We see the degree 6 and 12 polynomial lower approximations in Figure g They 
are both computed in less than 2 seconds. 


Example 12. As in Examples g and g consider 

p ■ s4 p(q, s) = s 3 + |s 2 + g 2 s + (g - \)q{q + |). 

The abscissa a p / of p' is not differentiable in two points, hence it is not a polynomial 
and it cannot be described perfectly by v d ' f° r finite d!. Let us choose d' — 8 and d — 6 
resp. d = 12. We observe in Figure [10] that Wq resp. w 12 is not everywhere a valid lower 
bound. Indeed, the set V = {q G Q \ a p (q) = a p '{q)} contains three points and for two 
of these (near q = —0.5 and q — 0), the approximation vg is not tight enough to ensure 
nQ(Z rt s) = Q. Consequently, Assumption g is violated. 
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Figure 10: Abscissa (black) and its polynomial lower approximations of degree 6 (gray, 
left) and 12 (gray, right) for Example 12 We observe that the approximations are not 


valid near q = —0.5 and q = 0, as Assumption [2] is violated. 


Example 13. In order to discuss another example for which V is a non-empty interval, 
consider the polynomial 


p : s p(q, s) = s 2 + (20g 2 — l)s + q + 


Here a p >{q) = 
particular v 2 


—10 q 2 + \ is a quadratic polynomial. Thus, Assumption [2l 
= dpi, and the lower approximations are valid, see Figure 11 


is fulfilled, in 


Example 14. As in Examples [4] and [9] consider the polynomial 
p: p(q, s) = s 3 + (qi + §)s 2 + q\s + q x q 2 . 


we see the outer approximations of degree d — 6 resp. d — 12 of the stabilizability region 
obtained for the choice d' = 8. A careful examination reveals that Assumption [ 2 ] is slightly 
violated here, yet this has no effect on the validity of the zero sublevel set approximation. 
Computing the degree 12 approximation takes a few minutes. 


We have Z rd t := {(q,x,y) G Z : x — Vd'{q) > 0} with Z given in Example W] In Figure 12 


5 Conclusion 


In this paper we continued our long haul research programme consisting of developing and 
applying semidehnite programming hierarchies for approximating potentially complicated 
objects (arising in optimization and control) with simple objects, namely polynomials of 
given degrees. The complicated object of interest here was the polynomial abscissa, which 
has low regularity, while being ubiquitous in linear systems control. 

In section [3] we described how to construct polynomial upper approximations to the abscissa 
with guarantees of L 1 convergence (or equivalently almost uniform convergence) on compact 
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Figure 11: Abscissa (black) and its polynomial lower approximations of degree 6 (gray, 
left) and 12 (gray, right) for Example 13 


sets. Constructing polynomial lower approximations with similar convergence guarantees 


has proved to be much more challenging. We proposed a first approach in Section 4.1 


using elementary symmetric functions which is quite general but also computationally 
challenging due to the introduction of many lifting variables. This motivated the study of a 
second approach in Section 442 using the GauB-Lucas theorem which is less computationally 
demanding, but unfortunately much more involved and subject to working assumptions. 


An interesting question that would deserve careful investigation is whether our L 1 con¬ 
vergence guarantees can be strengthened to L°°, i.e. to uniform convergence, since we 
know that the polynomial abscissa is continuous, and hence that it can be uniformly 
approximated by polynomials on compact sets. For this the semidehnit-e programming 
hierarchy should be modified accordingly. 
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